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We present our updated results of the nucleon-nucleon potential in quenched lattice QCD simula- 
tions with the plaquette gauge action and the Wilson quark action on the 32 4 (~ (4.4fm) 4 ) lattice. 
From the equal-time Bethe-Salpeter (BS) wave function, the NN potential is constructed through 
the Schrodinger-type equation. Resulting NN potential has all the qualitative features which phe- 
nomenological potentials commonly have: the repulsive core at short distance and the attractive 
well at medium and long distances. In the L — > °° limit, our NN potential is guaranteed to re- 
produce the scattering length obtained from the Liischer's formula. The quark mass dependence 
of the NN potential is studied with m K ~ 380,529,731 MeV. The results suggest that both the 
repulsive core at short distance and the attractive well at medium distance are enhanced in the 
light quark mass region. 
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1. Introduction 

The nuclear force is the essential ingredient in nuclear physics. Since H. Yukawa has intro- 
duced the pion 72 yeas ago |jlj], enormous efforts have been devoted to understand the origin of the 
nuclear force [gj. At medium to long distances (r^l.2 fm), the nuclear force is attractive, which 
is essential for the existence of bound nuclei. Furthermore, the nuclear force exhibits the repulsive 
core at short distance (r^0.7 fm), which is intimately related to the stability of heavy nuclei and 
neutron stars. The field theoretical treatment on the basis of the pion and heavy meson exchanges 
has reasonable success at long and medium distances. However, the repulsive core at short dis- 
tance, which was first introduced by R. Jastrow 56 years ago to explain the high energy behavior of 
the NN phase shifts [Q], has not been understood with firm theoretical ground. Since the nucleons 
begin to overlap at short distances, the quark and gluon structure of the nucleon is expected to play 
an important role. In this respect, the first principle lattice QCD calculation is desired to serve as 
the most powerful tool to attack the problems of the nuclear force. There exits a previous attempt 
to apply the formalism of the static qq potential to NN potential [Q]. However, the method which 
employs the static quarks is not faithful to the scattering data of nucleons. 

We have recently generalized the method developed by CP-PACS collaboration [§] for the 
nn scattering length and have carried out a first calculation of the NN potential from lattice QCD 

In our method, the Bethe-Salpeter wave function, whose asymptotic form leads to correct 
NN phase shift in the asymptotic regime, is introduced to construct the NN potential through the 
Schrodinger-type equation. The resulting NN potential has all the qualitative features required 
by phenomenology, i.e., the repulsive core at short distance and the attractive well at medium to 
long distances. In the report, we extend our previous calculations by increasing statistics and by 
introducing different quark masses corresponding to the pion masses, m % ~ 380,529,731 MeV. 
The results show that both the repulsive core and attractive well are enhanced in the light quark 
mass region. 

2. The formalism 

We begin with the Schrodinger-type equation, which is satisfied by the Bethe-Salpeter wave 
function for the NN system, 

(v 2 -* 2 W £ (r) = m N ^U(r,f)y E (r') (2.1) 

V E , a p(7) = lim (0\T [ Pa (x,t)n p (y,0)] \NN(E)) (2.2) 

where r = x — y. p a {x) and np(y) denote the standard local interpolating fields for nucleons, 

p a (x) = Babe (u T a (x)Cy 5 d b {x)) u CM {x), np(y) = e abc (u T a (y)Cy 5 d h {y)) d c> p(y), (2.3) 

where a,b and c denote the color indices, a and j8 denote the Dirac indices. C denotes the charge 
conjugation matrix. (For derivation of Eq. ( |2. 1[ ), see Refs. [j5L [7L H].) U(r,r'), which does not 
depend on E, plays a role of the non-local interaction kernel [Q]. The most general (off-shell) form 
of NN potential after imposing constraints arising from various symmetries is analyzed in Ref. |^]. 
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By applying the derivative expansion up to 0(V 2 ), we obtain 



U{r/) = [/ 1 (r,r / ) + (f 1 -f 2 )[/ TT (r,r') 

= ^('=%('=°)(r,rO + ^ (/ =% (/=1) (^) 
U {1 \r,r>) = V {1) -8{x-x') 



C + (cti • o 2 )VX> + S 12 V^ + (L ■ S)V^ + L, a 2 -l} + {a! • V, a 2 • v} r , ip 



(2.4) 



* V (/) (r) + (dj • dz)V*> (r) + S 12 ^» + (I • (r) + 0{ V 



,0 



(2.5) 



where ! ^ | ( / = ) = (1 - ?i • T 2 )/4 and J 2 ^ 1 ) = (3 + T\ ■ f 2 )/4 denote the projection matrices to / = 
and 1=1 subspaces, respectively. {*, *} denotes the anti-commutator. Sn = (oi • r)(a 2 • r)/r 2 — 
dy -d 2 X = -if x V and S = (di + a 2 ) /2. ^ , , , v£ } and V$ are functions of r 2 , V 2 , and 
L 2 . We combine the 1st and the 2nd terms in the last line as V^p(r) = Vq (r) + (ai • CT 2 )Va (r), 
and refer to vi ! \r) as the central force. Vj\r) and Vyi (r) are refereed to as the tensor force and 
the LS force, respectively. These three forces play major roles in conventional nuclear physics. 
In QCD, the closest concept to the quantum mechanical wave function is provided by the 

represents a probability 



equal-time Bethe-Salpeter (BS) wave function Eq. (|2.2|). Note that Eq. (|2 
amplitude to find three quarks at x and another three quarks at y. It is possible to show that Eq. (| 
has a proper asymptotic behavior at \x — y \ — > oo [M]. For example, in the channel, we have 



JSo(k) 



sin (kr+8o(k)) 
kr 



(2.6) 



The BS wave function is obtained from the large t behavior of the "four-point" correlator of the 
nucleon 



(0\T[p(x,t)n(x,t)W(t = t )]\0) = Y,(0\p{x,0)n(y,0)\n)e 



-iE n (t-t ) 



(n|W(t = 0)|0). (2.7) 



Here, W(t) = P(t)N(t) represents the wall source, where P(t) and N(t) are defined as Eq. (23) with 
the quark fields q(x,t) replaced by Q(t) = L^^(x,f). Since contributions from all excited states are 
exponentially suppressed in the large t region, we are left with the BS wave function for the ground 
state. In this report, we consider only 1 Sq and 3 Si channels whose BS wave functions are 



l 



r3 

ReO 1 - x 



ZH ReO u x 



p a (R-r+X)n p (X) NN(E) 



Pa (R-r + X)n p (X) NN(E)), 



(2.8) 



where the summations over R G O are performed for cubic transformation group. The summations 
for X are performed to select the zero total spatial momenta. 

First we consider the Schrodinger equation in l So channel. Owing to the identical nature of the 
nucleons, two nucleon system in l So channel is iso- vector. Since the contributions from the tensor 
force and the LS force vanish in this channel, we are left with the following Schrodinger equation, 



1 rt 
2/1 



V 2 + V^- l \r)^(r; l S ) = ^(r; 1 



So), 



(2.9) 



3 



Lattice QCD approach to nuclear force 



N. Ishii 



K 


Aconf 


m n [MeV] 


ra N [MeV] 


t-to ECS )[MeV] 


£( 3 5i) [MeV] 


0.1640 
0.1665 
0.1678 


1000 
2000 
2021 


732.1(4) 
529.0(4) 
379.7(9) 


1558.4(63) 
1333.8(82) 
1196.6(32) 


7 -0.400(83) 
6 -0.509(94) 
5 -0.675(264) 


-0.480(97) 

-0.560(114) 

-0.968(374) 



Table 1: The number of gauge configurations N CO nf, the pion mass ?n n , the nucleon mass m^, time-slice 
t — to on which BS wave functions are measured, and the non-relativistic energies E = k 2 /(2n) for % and 
3 Si channels. 



where pi = m^/2 denotes the reduced mass of the nucleon. k plays the role of the "asymptotic 
momentum". Since V^ =l) {r) = V (/=1) (r) - 3V^ I=l) (r) is an ordinary function, which does not 
involve a derivative nor matrix structure, we arrange the Schrodinger equation to obtain 



c W 2^ 2^ ^(r; 1 ^) 

Next, we consider the Schrodinger equation in 3 Si channel. In this case, two nucleon system 
is iso-scalar. Unlike 1 Sq case, this channel receives a non- vanishing contribution from the tensor 
force, which provides a coupling to 3 D\ channel. 3 D\ channel receives a contribution from the LS 
force. In this way, we have three unknowns, i.e., VQ~°\r), Vj~°\r) and V^~°\r) with two equa- 
tions, i.e., 3 Si and 3 D X . (V^ /=0) (r) = V (/=0) (r) + vi /=0) (r) should not be confused with V<f =1) (r).) 
To obtain these three forces exactly, we need one more equation such as the Schrodinger equation 
in 3 D2 channel. In this report, we do not pursue this direction. Instead, we adopt the same pro- 
cedure as the 1 Sq channel. This leads to the so-called "effective central force" V£ (r) which takes 
into account the 3 D\ channel indirectly through the tensor force. 



3. The lattice QCD result 

We employ the standard plaquette action on 32 4 lattice with /3 = 5.7 to generate quenched 
gauge configurations. The gauge configurations are picked up every 200 sweeps after skipping 
3000 sweeps for thermalization. Quark propagators are generated by employing the standard 
Wilson quark action with K = 0.1640,0.1665 and 0.1678. The scale unit l/a = 1.44(2) GeV 



(a ~ 0.137 fm) is introduced from the rho meson mass in the chiral limit QlOp . The physical size 
of our lattice corresponds to L ~ 4.4 fm. The number of gauge configurations A^ CO nf, the pion mass 
m K , the nucleon mass m^ are summarized in Table [[]. (For K = 0.1678, 24 gauge configurations 
are identified as exceptional configurations, which are not used in the calculations.) 

The periodic (Dirichlet) boundary condition is imposed on the quark fields along the spatial 
(temporal) direction. We adopt the wall source on the time-slice t = to = 5. The BS wave functions 
are measured on the time-slice t — to = 7,6,5 for K = 0.1640,0.1646,0.1678, respectively. The 
ground state saturation is examined by t-dependence of the NN potential. We employ the nearest 
neighbor representation of the discretized Laplacian as V 2 f(x) = £ 3 =1 {f(x + £?£,•) +f(x — aef)} — 
6/ (3c), where e, denotes the unit vector along the i-th coordinate axis. BS wave functions are fully 
measured for |r|^$0.7 fm, where rapid changes of BS wave function and NN potential are expected. 
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Figure 1: NN wave functions in 'So and Si channels (left), and NN potentials (right), i.e., central force in 



'So channel and effective central force in 3 Si channel for K ■ 
of the wave function <j)(x,y,z = 0;' So). 



0.1665. The inset of the left figure is a 3D plot 



Since the changes are rather modest for |r|^0.7 fm, the measurement of BS wave functions is re- 
stricted on the coordinate axes and their nearest neighbors to reduce the calculational cost. The 
"asymptotic momentum" k 2 is obtained by fitting the BS wave function: We use the Green's func- 
tion of the Helmholtz equation, 



G{r;k 2 



1 r- 1 exp(/27T«-r/L) 
tfj- (2n/L) 2 n 2 -B- 



(3.1) 



as the fit function by regarding the overall numerical factor and k 2 as fit parameters. (An appropriate 
regularization is assumed in this representation of Green's function.) The fits are performed outside 
of the range of NN interaction, which is determined by examining V 2v F(*) /^(x) [§]. 

Fig. [T] (left) shows the BS wave functions in l So and 3 Si channels for K = 0.1665. The sup- 
pression of the wave function in the region r<0.5 fm indicate the existence of repulsion at short 
distance. Fig. [I] (right) shows the reconstructed NN potentials for K = 0.1665, i.e., the central force 
for 'So channel and the effective central force for 3 Si channel. (See Table ffl, for the values of the 
non-relativistic energies E = k 2 / in Eq. (2. 1C ).) We see that our NN potentials have repulsive 
cores of 500 — 600 MeV in the short distance region (r<0.5 fm) and attractions of about 30 MeV 
in the medium distance region (0.5^jr<1.2 fm). Both the repulsive core at short distance and the 
attractive well at medium distance are weaker than those expected phenomenologically. This is 
due to the heavy quark mass in our simulations In the light quark mass region, pion can propagate 
longer distance, which is expected to enhance the attraction at medium and long distance. 

To see the quark mass dependence in a quantitative manner, we show the NN potentials for 
three different quark masses in 'So channel in Fig. ||. As the quark mass decreases, the repulsive 
core at short distance is enhanced rapidly, whereas the attraction at medium distance is modestly 
enhanced. This indicates that it is important to perform the lattice QCD calculation in the lighter 
quark mass region in order to compare our result with experimental data. 

Several comments are in order here, 
(i) The net interaction is attractive even in the presence of the repulsive core. Indeed, Luscher's 
finite volume method leads to the attractive scattering length, i.e., ao('So) = 0.115(26), 0.126(25), 
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Figure 2: Central forces in So channel for three quark masses. 



0.159(66) fmand a ( 3 Si) = 0.140(31), 0.140(31), 0.252(104) fm for k = 0.1640,0.1665,0.1678, 
respectively. The attractive nature of our potential is qualitatively understood by the Born approxi- 
mation formula for the scattering length ciq ~ — § Vc{r)r 2 dr. Owing to the volume factor r 2 dr, 
the attraction at medium distance overcomes the repulsive core at short distance, 
(ii) There is a considerable discrepancy between the above scattering lengths and the empirical val- 



20 fm and ai ( 3 5i) ~ —5 fm. This is attributed to the heavy quark mass 



ues, i.e., c?| ) exp) ( 1 5 , ) 

employed in our simulations. If we can get closer to the physical quark mass, there appears an 
"unitary region" where the NN scattering length becomes singular as a function of the quark mass 
and changes sign JTT| , O ]. The singular point is related to the threshold of bound state formation. 
This is why the physical scattering length is positively large in the 1 Sq channel (no bound state) and 
is negatively large in the 3 S\ channel (deuteron bound state) . 



4. Summary 

We have extended our previous results of the nuclear force on the lattice by increasing statistics 
and adopting different quark masses. The NN potentials in the l So and 3 Si channels have all the 
qualitative features which phenomenological NN potentials commonly have, i.e., the repulsive core 
at short distance and attractive well at medium to long distances. The quark mass dependence of the 
NN potential shows that the repulsive core at short distance is enhanced rapidly, and the attraction at 
medium distance is modestly enhanced. These results suggest that, in order to compare our results 
with the experimental data, it is important to perform the lattice QCD Monte Carlo calculation in 
the lighter quark mass region. 
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Although the BS wave functions are proved to have the universal behavior as Eq. (2.6) at 
large distance, its short distance behavior is afflicted with the operator dependence. We can avoid 
this subtlety by resorting to the inverse scattering theory, which guarantees the unique existence 
of energy-independent, local, hermitian potential. Studies along this line together with the NN 
potential measured in different E and the derivation of the tensor force will be reported elsewhere 
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